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ABSTRACT 

We report on the results from Suzaku broadband X-ray observations of the galactic binary source LS 5039. 
The Suzaku data, which have continuous coverage of more than one orbital period, show strong modulation of 
the X-ray emission at the orbital period of this TeV gamma-ray emitting system. The X-ray emission shows 
a minimum at orbital phase ^0.1, close to the so-called superior conjunction of the compact object, and a 
maximum at phase ~ 0.7, very close to the inferior conjunction of the compact object. The X-ray spectral 
data up to 70 keV are described by a hard power-law with a phase-dependent photon index which varies 
within r ~ 1.45- 1.61. The amplitude of the flux variation is a factor of 2.5, but is significantly less than 
that of the factor ^8 variation in the TeV flux. Otherwise the two light curves are similar, but not identical. 
Although periodic X-ray emission has been found from many galactic binary systems, the Suzaku result implies 
a phenomenon different from the "standard" origin of X-rays related to the emission of the hot accretion plasma 
formed around the compact companion object. The X-ray radiation of LS 5039 is likely to be linked to very- 
high-energy electrons which are also responsible for the TeV gamma-ray emission. While the gamma-rays 
are the result of inverse Compton scattering by electrons on optical stellar photons, X-rays are produced via 
synchrotron radiation. Yet, while the modulation of the TeV gamma-ray signal can be naturally explained 
by the photon-photon pair production and anisotropic inverse Compton scattering, the observed modulation 
of synchrotron X-rays requires an additional process, the most natural one being adiabatic expansion in the 
radiation production region. 

Subject headings: acceleration of particles — X-rays: individual (LS 5039) — X-rays: binaries 
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1. INTRODUCTION 



LS 5039 is a high-mass X-ra v binary (iMotch et al l 1 19971) 
with extended radio emission dParedes et alJ 120001 12002). 
This system is formed by a main sequence O type star and 
a compact object of dispu ted nature that has b een claimed to 
be both a black hole (e.g..lCasares et al. 2005) and a neutron 
star/pulsar (e.g.. iMartocchiaet al. 2005: iDubusI 120061) . The 
compact object is moving around the companion star in a 
moderately elliptic orbit (ecce ntricity e = 0.35) w ith an orbital 
period of f orb = 3.90 60 days dCasare s et al. 2005). 

As summarized in Bosch-Ramon et all (12007b . LS 5039 has 
been observed several times in the X-ray energy band for 
limited phases in the orbital period. Flux variations on time 
scale of days and sometimes on much shorter timescales have 
been reported. The spectrum was always well represented 
by a power-law model with a photon index ranging T = 1 .4- 
1.6 up to ~ 10 keV, with fluxes changing moderately around 
~ 1 x 10 -11 erg cm" 2 s" 1 . Softer spectra and larger fluxes had 
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been also inferred from RXTE observations, although back- 
grou nd contamination was prob ably behind these differences 
(see iBosch-Ramon et al.1 120051) . Also, Chandra data taken 
in 2 002 and 2005 showed sp ectra significantly harder than 
1 .5 dBosch-Ramon et aT]|2007l) . but such a hard spectrum was 
probably an artifact produced by photon pile-up. Recently, 
iHoffmann etail (120091) reported the results of INTEGRAL ob- 
servations in hard X-rays. The source was detected at en- 
ergies between 25 and 60 keV. The source was detected at 
energies between 25 and 60 keV. The flux was estimated to 
be (3.54 ±2.30) x 10" 11 erg cm" 2 s" 1 (90 % confidence level) 
around the inferior conjunction (INFC) of the compact object, 
and a flux upper limit of 1 .45 x 10 -11 erg cm" 2 s" 1 (90 % con- 
fidence level) was derived at the superior conjunction of the 
compact object (SUPC). 

LS 5039 has also been detect ed in very-high-energy ( VHE: 
E > 0.1 TeV) gamma-rays ( Aharoni an et alJ [2005al) . ex- 
hibiting a periodic si gnal modulated with the orbital period 
(Aharonian et alj|2006l) . There are two other binary systems 
with robust detections i n the VHE range: PS R B 1259-63 
dAharonian et alJ l2005bl) and LS I +61°303 dAlbert et alJ 
2006; lAcciari et alJ 12008). Evidence for TeV emissio n 
has been found also in Cygnus X-l (Alb ert et all l200l . 
PSR B 1259-63 is a clear case of a high-mass binary sys - 
tem containing a non-accreting pulsar (Johnston et al. 1992), 
whereas Cygnus X- 1 is a well known accreting black hole 
system dBoltonlll972l) . The nature of the compact object in 
LS 5039 is not yet established, and the origin of the VHE 
emitting electrons is unclear. They may be related to a pulsar 
wind or to a black hole with a (sub)relativistic jet. In the stan- 
dard pulsar-wind scenario, the severe photon-photon absorp- 
tion makes the explanation of the detection of the VHE radia- 
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tion problematic, at least at the position corresponding to the 
orbital phase <f> ~ 0.0 (see Figs. 16 and 4 from Sierpowska- 
Bartosik & Torres 2008 and Dubus et al. 2008 and compare 
with Fig. 5 in Aharonian et al. 2006), in which the emitter 
is expected to be located between the compact object and the 
star. However, particles may be accelerated in a relativistic 
outflo w formed at the interact ion of the pulsar and the stellar 
winds (Bogoval ov et alj |2008) and radiate far from the com- 
pact object, making the pulsar wind scenario a viable option. 
In the microquasar scenario, the lack of accretion features in 
the X-ray spectrum may be a problem unless the bulk of the 
accretion power is released in the form of kinetic energy of 
the outflow, rather than thermal e mission during acc retion, 
as in the case of SS 433 (see e.g., iMarshall et alj|2002l) . At 
this stage, we cannot give a preference to any of these scenar- 
ios, but new data, in particular those obtained with the Suzaku 
satellite, allow us to make an important step towards the un- 
derstanding of the nature of the non-thermal processes of ac- 
celeration and radiation in this mysterious object. 

2. OBSERVATION 

The temporal and spectral characteristics of the X-ray emis- 
sion from LS 5039 along the orbit should provide impor- 
tant clues for understanding the acceleration/radiation pro- 
cesses in this source. The fact that all previous X-ray ob- 
servations of this object have incomplete coverage of the 
orbital period, or suffered from background contamination, 
is therefore rather unsatisfactory. This motivated our long, 
~ 200 ks observatio n with the Suzaku X-ray observatory 
(Mitsu da et al.ll2007h . which gives us unprecedented cover- 
age of more than one orbital period, continuously from 2007 
September 9 to 1 5 (see T able Q]). Suza ku has four sets of 
X-ray telescopes (Serlemitsos et al. 2007) each with a focal- 
plane X-ray CCD ca mera (X-ray Imaging Spectrometer(XIS); 
Kovar ria et al.l l2007) that are sensitive in the energy range of 
0.3-12 keV. Three of the XIS detectors (XIS0, 2 and 3) have 
front-illuminated (FI) CCDs, whereas XIS1 utilizes a back- 
illuminated (BI) CCD. The merit of the BI CCD is its im- 
proved sensitivity in the soft X-ray energy band below 1 keV. 
Suzaku contains also a non-ima ging c ollimated Hard X-ray 
Detector (HXD; Takahas hTeUd]l2007l Ifokubun et al.ll2007h . 
which covers the 10-600 keV energy band with Si PIN pho- 
todiodes (10-70 keV) and GSO scintillation detectors (40- 
600 keV). Suzaku has two default pointing positions, the XIS 
nominal position and the HXD nominal position. In this ob- 
servation, we used the HXD nominal position, in which the ef- 
fective area of HXD is maximized, whereas that of the XIS is 
reduced to on average ~ 88%. Results from XIS2 are not re- 
ported here since it has not been in operation since an anomaly 
in November 2006. In addition, we do not describe in detail 
the analysis of HXD-GSO data, since the HXD-GSO detected 
no significant signal from the source. 

3. DATA REDUCTION 

We used data sets processed using the software of the 
Suzaku data processing pipeline (version 2.1.6.16). Reduc- 
tion and analysis of the data were performed following the 
standard procedure using the HEADAS v6.4 software pack- 
age, and spectral fitting was performed with XSPEC v. 1 1.3.2. 

For the XIS data analysis, we accumulated cleaned events 
over good time intervals that were selected by removing 
spacecraft passages through the South Atlantic Anomaly 
(SAA). Further, we screened the data with the following crite- 
ria — (1) cut-off rigidity is larger than 6 GV and (2) elevation 
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FIG. 1. — The HXD-PIN spectrum. The black, red, blue, green, and ma- 
genta points show the count rates based on the raw data, from the modeling 
of NXB, the CXB component, the GRXE component, and the background- 
subtracted data, respectively. The background here includes contributions 
from the NXB, CXB and GRXE (see §2 for details). 

angle from the Earth's rim is larger than 5°. The source pho- 
tons were accumulated from a circular region with a radius of 
3'. The background region was chosen in the same field of 
view with the same radius and an offset of 9' from the source 
region. We have co-added the data from the two FI-CCDs 
(XIS0 and XIS3) to increase statistics. The response (RMF) 
files and the auxiliary response (ARF) files used in this paper 
were produced using xisrmf gen and xissimarf gen, re- 
spectively. 

For the HXD data,"uncleaned event files" were screened 
with the standard event screening criteria: the cut-off rigid- 
ity is larger than 6 GV, the elapsed time after the pas- 
sage of the SAA is more than 500 s and the time to 
the next SAA passage is more than 180 s, high volt- 
ages from all eight HV units are within the normal range 
and the elevation angle from the Earth's rim is more 
than 5°. We also discarded telemetry-saturated time inter- 
vals. In the spectral analysis in §3, we used the response 
file for a point-like source at the HXD-nominal position 
(ae_hxd_pinhxnome4_20070914 . rsp), which is re- 
leased as a part of CALDB (Suzaku calibration data base). 

The HXD-PIN spectrum is dominated by the time-variable 
instrumental background ( non-X-ray background (NXB) ) 
induced by cosmic rays and trapped charged particles in the 
satellite orbit. To estimate the instrumental background com- 
ponent, we used the the time-dependent NXB event files re- 
leased by th e HXD instrument team , whose reproducibility is 
reported by Fukaza wa et alj (12009!) . In order to estimate the 
systematic uncertainty of the NXB model during our obser- 
vation, we compare the NXB model spectrum during Earth 
occultation with the observed spectrum of the same time in- 
terval. The event selection criteria for this study are the same 
as those of the cleaned event except for the criterion on Earth 
elevation angle, which was chosen to to be less than -5°. The 
estimated uncertainty obta ined is ~ 3%, which is consistent 
with the values reported by Fukazawa et al. (20Q9|). 

Another component of the HXD-PIN background is the 
cosmic X-ray background (CXB). In our analysis, we as- 
sumed the CXB spectrum reported by dGruber et alj fl999): 

7(e) = 7.9 Ekev 9 exp ) P h s "' keV_1 cm ~ 2 str "' » (!) 

where e^v = e/keV and e p = 41.1. The CXB spec- 
trum observed with HXD-PIN was simulated by us- 
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TABLE 1 
Log of Suzaku Observations 



Obs. ID 


Coord. (J2000) 


Exposure 


Date 




RA, DEC 


XIS/HXD 




402015010 


18 h 26 ra 15 s .l,-14 d 50 ra 54 s .2 


203 ks/181 ks 


09-15/09/2007 



Note. — The exposure time is the net integration time after standard data screening for the XIS and HXD-PIN. 

ing a PIN response file for isotropic diffuse emission 2 
(ae_hxd_pinflate4_2 007 914 . rsp) and added to 
the NXB spectrum. Based on this approach, the contribution 
from the CXB flux is ~ 5% of the NXB. 

Since LS 5039 is located close to the Galactic plane, the 
contribution from the Galactic ridge X-ray emission (GRXE) 
must be examined, especially for HXD-PIN spectra. In order 
to model the shape of the GRXE, data from Suzaku obser- 
vations of the Galactic ridge region (ObsID: 500009010 and 
500009020) were analyzed. The Suzaku spectrum from 3 keV 
to 50 keV can be well fitted with the Raymond-Smith plasma 
with a temperature of kT = 2.2 ± 0.8 keV and a power-law 
function with T = L92^;f . Although we also tried a power 
law with expone ntial cutoff, following the results from the 
INTEGRAL IBIS dKrivonos et alJl2007l) . it turned out that the 
assumption on the spectral shape of the GRXE has negligible 
effect on the spectral parameters of LS 5039. The normaliza- 
tion of the GRXE spectrum component in the HXD-PIN spec- 
trum of LS 5039 is determined from the XIS spectrum of the 
LS 5039 observation by excluding an encircled region with a 
radius of 4.5' centered on the LS 5039 location. The flux of 
the GRXE is estimated to be ~ 40% of the contribution from 
the CXB. In Figure Q] we show the time averaged HXD-PIN 
spectrum plotted together with models for the NXB, CXB and 
GRXE. 




4. ANALYSIS AND RESULTS 
4.1. Temporal analysis 

The light curve obtained from the XIS detector is shown in 
the top panel of Figure[2] The continuous coverage in X-rays, 
longer than the orbital period of the LS 5039 system, reveals 
a smooth variation of a factor 2 in the 1-10 keV count rate. 
The light curve is drawn over two orbital periods. The orbital 
phase is calculated with the period of 3.90603 days, and <f> = 
with reference epoch T n ( HJD- 2400000.5 = 51942.59) taken 
from|Casares et al. (2005). The light curve from phase 0= 1.0 
to 1.5, which was obtained in the last part of the observation, 
smoothly overlaps with the one obtained at the beginning of 
the observation (0 = 0.0-0.5). 

In the middle panel of Figure [2] we present the light curve 
obtained with the HXD-PIN for the energy range 15^40 keV. 
Although the statistical errors are larger, the modulation be- 
havior is similar to that of the XIS. The amplitude of the mod- 
ulation is roughly the same between the XIS and HXD-PIN, 
indicating small changes of spectral shape depending on or- 
bital phase. The spectral parameters obtained for each orbital 
phase are reported in the following section. 

The light curves obtained with Suzaku show that the X-ray 
flux minimum appears around phase 0.0-0.3 and it reaches 
maximum around phase 0.5-0.8. In order to quantify the am- 
plitude of the flux variations, we fitted the XIS light curve 
with a simple sinusoidal function. Due to structures in the 
light curve, the fit converges with large chi-square (x^(^) = 
4.92 (121)). However, the general trend is well represented 
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Orbital Phase [0] 

FIG. 2. — Orbital light curves observed for LS 5039. (Top) Suzaku XIS 
count rate in the 1-10 keV band continuously obtained for a 200 ks duration 
(filled circles), which covers one orbit and a half. Overlaid is the same light 
curve but shifted by one orbital period (open circles). (Middle) Suzaku PIN 
count rate in the 15^10 keV band. Background components consisting of 
NXB, CXB and GRXE are subtracted. (Bottom) the light curve of integral 
fluxes at energies E> 1 TeV vs. orbital phase (<f>) obtained on run-by-ru n 
basis from the HESS data (2004 to 2005) reported by Aharonia netal] f2006l) . 



by a sinusoidal function: 



I(4>) = 0.24 sin(</>- 0. 1 3) + 0.64 counts s" 



(2) 



where I(<fi) is the count rate as a function of phase <\>. The 
ratios between the minimum and maximum count rates are 
2.21 counts s" 1 for XIS and 2.02 ^ f 9 counts s" 1 for 
HXD-PIN. Structures of the X-ray and hard X-ray light curves 
are similar to that discovered in the phase diagram of inte- 
gral fluxes at energies > 1 TeV obtai ned on a run-by-run ba - 
sis from HESS data (2004 to 2005) (lAharonian e t al. 2006). 
The temporal X-ray behavior was already suggested by RXTE 
data, as well as by a compilation of all the previous X-ray 
data obtained with imagin g instruments (Bosch- Ramon et al] 
120051; IZabalza et al.ll2008l) . 

In addition to the continuously changing component with 
respect to the orbital phase, short timescale structures are 
found around <fi = 0.48 and cf) = 0.7. The unabsorbed flux 
changes about a 30% in A<j> ~ 0.05 (= 4.7 hour). A significant 
dip can be seen around <fi = 1 . 35 in the top panel of Figure|2] In 
comparison with the data at around </> = 0.35 obtained in the 
first half of the observation, the flux decreased ^50% only 
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FIG. 3. — Orbital light curves of the photon indices, Nn, and unabsorbed 
fluxes in the energy range of 1-10 keV obtained by fitting with an absorbed 
power-law model. The phase interval width is A<j> = 0.1. Two full phase pe- 
riods are shown for clarity. The blue solid lines show periastron and apastron 
phase and red dashed lines show SUPC and INFC of the compact object. 

in this phase. The time duration of this dip corresponds to 
Atfi ~ 0.03. These structures may reflect features of the (pos- 
sibly changing) environment of the X-ray emitting region, it 
is therefore of importance to test with further observations if 
these are persistent features. 

4.2. Spectral Analysis 

Firstly we study time-resolved (phase-resolved) X-ray 
spectra. The data are divided into data segments with respect 
to the assigned phase, and model fitting is performed for XIS 
spectra for each segment with = 0.1. 

A single power-law function with photoelectric absorption, 
provides a good fit for all the segments. In order to study the 
possible changes of the amount of photoelectric absorption, 
we here fit the data with A^h free. The best-fit parameters are 
presented in Table [2] The derived values of the photon in- 
dex and absorpti on column density are consistent with previ- 
ous o bservations (Martoc chia et al.l2005tlBosch-Ramon et al.l 
120071) . When we fix the A^h to the value obtained from the time 
averaged spectrum, resultant photon indices stay same within 
statistical error. 

The photon index (T) values are plotted as a function of 
orbital phase in the top panel of Figure [3] The spectral 
shape varied such that the spectrum is steep around SUPC 
(r ~ 1 .61) and becomes hard {V ~ 1 .45) around apastron. The 
modulation behavior of T is somewhat different from that ob- 
served using HESS in the VHE range. The amplitude of the 
variation is ±0.1, which i s much smaller than th e change of 
±0.6 in the VHE region (lAharonian et al.|[2006b . The 1-10 

'H<t> = 
(0 = 0.6-0.7). 

In all the data segments, the source is significantly detected 
with the HXD-PIN, indicating that hard X-ray emission ex- 
tends at least up to 70 keV. Note also that although the XIS and 
HXD-PIN spectra do not overlap, they seem to be smoothly 
connected in the gap between 10 and 15 keV. 

To study the shape of the spectrum above 10 keV, the XIS 
spectra and the PIN spectrum in the range 15-70 keV after 



keV flux changes from (5.18 ±0.03) x 10 12 erg cm 2 s 
0.1-0.2) to (12.05 ±0.02) x 10" I2 erg cm" 2 s"' 



subtraction of background (NXB + CXB + GRXE) are jointly 
fitted (Figure [4] Top). The time-averaged spectra are well 
represented by an absorbed power-law model with T= 1.51 
±0.02 with reduced xl = 0.99 (235 de grees of freedom). We 
find no cutoff structure in the energy range of the HXD-PIN. 
The spectra within the phase intervals [0.616 < <f> < 0.816] 
and [(f> < 0.158 & 0.958 < </)], which correspond to the INFC 
and SUPC, respectively, are also shown in Figure[4] The best- 
fit parameters are presented in Table [3] 
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FIG. 4. — Top: The time-averaged Suzaku (XIS+HXD) spectrum (Top) and 
spectra obtained in INFC (Middle) and SUPC (Bottom). 



Althought earlier observations by RXTE suggested the pres- 
ence of an iron emission line at 6.7 keV (Ribo et al. 1999), 
later observation s by Chandra and XMM c ould not find evi- 
dence of it (e.g. Bosch- Ramon et al]|2005l) . A careful study 
of new and longer RXTE observations, using slew data to ac- 
count for background emission, revealed that the earlierly re- 
ported 6.7 keV emission line is likely a background feature 
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TABLE 2 

Results of Suzaku XIS spectral fitting 



Orbital Phase Intervals 


Photon index 


^Vh (wcibs) 


Flux [1-10 keV] 


Y 2 (U) 






10 21 cm" 2 


10 ^ 2 erg cm 2 s 


TNFC (0 45 < eh < 91 




7 82 + 1 8 


1 n 78 4- o 05 


0.97 (1 1 1) 


SIIPfYrji < 45 9 <r nil 


1 55+ 02 

' --■>- • ' \J .\Jj— 


7 77 _|_ o.20 


f. 79 + 04 


0.94 (1 1 1) 


All Phases (time averaged) 


1.51±0.01 


7.71±0.2 


8.07 ± 0.03 


0.81 (118) 


0.0-0.1 


1.57 ±0.04 


6.86 ± 0.44 


5.62 ± 0.03 


0.99 (111) 


0.1-0.2 


1.61 ±0.04 


7.67 ± 0.46 


5.18 ± 0.03 


1.04(111) 


0.2-0.3 


1.51 ±0.03 


7.23 ± 0.37 


5.67 ± 0.02 


1.13(107) 


0.3-0.4 


1.49 ±0.03 


7.42 ± 0.43 


7.34 ± 0.03 


1.01 (111) 


0.4-0.5 


1.45 ±0.02 


7.63 ± 0.26 


9.73 ± 0.01 


0.96(111) 


0.5-0.6 


1.46 ±0.03 


7.72 ± 0.45 


9.95 ± 0.02 


1.15 (105) 


0.6-0.7 


1.46 ±0.03 


7.82 ±0.37 


12.05 ± 0.02 


0.77 (118) 


0.7-0.8 


1.51 ±0.02 


7.61 ±0.33 


11.27 ±0.02 


1.10(108) 


0.8-0.9 


1.52 ±0.04 


7.63 ± 0.47 


10.29 ± 0.03 


0.87 (111) 


0.9-1.0 


1.59 ±0.03 


8.35 ±0.41 


7.84 ± 0.02 


1.06(111) 



Note. — Fitting Suzaku XIS0 +XIS3 spectrum of LS5039 by a power-law with photoelectric absorption in 0.6-10 keV. Photon index T, absorbing column density 
N}i, and the 1 — 10 keV flux F[_io (corrected for absorption) are shown with 68.3 % (1 cr) error. The orbital phase </> is calculated from the ephemeris of 
Casaresetal. (2005). 



TABLE 3 

Results of LS 5039 XIS+HXD-PIN spectral fitting 



Orbital Phase Intervals 


Photon index 


N H (wabs) 
Xl0 21 cm 2 


X>) 


All Phases (time averaged) 


1.51 ±0.02 


7.7 ±0.2 


0.95 (175 d.o.f.) 


INFC (0.616 <cf>< 0.816) 


1.49 ±0.04 


7.9 ± 0.4 


0.87 


SUPC (0 < 0.158,0.958 < <f> < 1.158) 


1.55 ±0.05 


6.9 ± 0.5 


1.06 



Note. — Power-law fitting to the XIS/HXD-PIN spectrum. The errors of spectral 

dBosch-Ramon et aLI 120051) . The Suzaku data confirm this 
result. In an attempt to find the possible signature of iron 
emission lines from LS 5039, we analyzed the phase-averaged 
spectrum. The upper limits on iron line structures are deter- 
mined by fitting a Gaussian at various energies and line widths 
at which Fe emission might be expected. The power-law con- 
tinuum model parameters are fixed with the best fit values and 
a Gaussian component is added to the power-law function. 
The central energy of the Gaussian line is swept from 6.0 keV 
to 7.1 keV in steps of 0.1 keV. Line widths are changed from 
0.01 keV to 0.09 keV in steps of 0.01 keV, together with lines 
with larger widths of 0.15 and 0.20 keV. An equivalent width 
is determined at each grid point. The resulting upper-limit on 
the equivalent width is 40 eV with 90% confidence level. 

5. DISCUSSION 

The X-ray emission observed with Suzaku is character- 
ized by (1) a hard power law with T ~ 1.5 extending from 
soft X-rays to ~ 70 keV, (2) clear orbital modulation in flux 
and photon index, (3) a moderate X-ray luminosity of Lx ~ 
10 33 (-r^) 2 erg s" 1 , (4) a small and constant absorbing column 
density, and (5) a lack of detectable emission lines. 

Although variable X-ray emission has been found from 
more than two hundred galactic binary systems, the Suzaku 
data hardly can be explained within the "standard accretion" 
scenario where X-rays are produced by a hot thermal (comp- 
tonized) accretion plasma around the compact object. The 
lack of X-ray emission lines (at the level of sensitivity of 
Suzaku) as well as the hard £ _1 5 type energy spectrum of the 
X-ray continuum, extending from soft X-rays up to 70 keV, 



parameters are at the 90% confidence level. 

favors a non-thermal origin of the X-rays. This conclusion is 
supported by the general similarities between the properties 
of the observed X-rays and TeV gamma-rays. Namely, both 
radiation components require a rather hard energy distribu- 
tion of parent electrons with a power-law index of » 2. This 
directly follows from the photon index of the synchrotron ra- 
diation r w 1.5, and agrees quite well with the currently most 
favored interpretation of the TeV gamma-rays, in which they 
would be produced by IC scattering off the anisotropic photon 
field of the massive companion star. 

Assuming that the TeV gamma-ray production region is lo- 
cated at a distance from the companion star of d ~ 2 x 10 12 cm 
(i.e. the binary system size), and taking into account that 
gamma-rays are produced in the deep Klein-Nishina (KN) 
regime with significantly suppressed cross-section, for the 
well known luminosity of the optical star L~7 x 10 38 ergs _1 , 
one can estimate quite robustly the strength of the magnetic 
field in the emission region. The numerical calculations show 
that the field should be around a few Gauss (see e.g. Fig. 
5). For such a magnetic field strength, the energy inter- 
vals of electrons responsible for the two emission compo- 
nents overlap substantially, as shown in Figure [5] There- 
fore, we are most likely dealing with the same population of 
parent electrons, which should be located at large distances 
from the compact object, in the system periphery, to pre- 
vent the severe absorption of the TeV radiation and the sub- 
sequent intense emission from the pair-created seco ndaries 
(Khan gulyaneta l. 2008a: iBosch-Ramon et aTll2008al) . 

It should be noted that the observed X-ray emission is very 
difficult to explain as synchrotron emission produced by sec- 
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FIG. 5. — The radiative cooling times as a function of electron energy. 
The blue line corresponds to IC losses at the distance d = 2 X 10 12 cm from 
the optical star (the star luminously and temperature were assumed to be 
L, = 7 X 10 38 erg/s and T t = 3.8 X 10 4 K, respectively). The red line corre- 
sponds to the synchrotron losses for B-field B = 3 G. The filled regions repro- 
duce the electron energy intervals relevant to the HESS (yellow), Fermi (light 
blue) and Suzaku (green) energy domains. Formally, two radiation chan- 
nels, synchrotron radiation of very high energy electrons (light green) and IC 
scattering of low energy electrons in the Thompson regime (dark green) can 
produce X-ray photons in the Suzaku energy domain. 

ondary (pair-produced) electron and positrons. Since the pair 
production cross-section has strong energy dependence with 
a distinct maximum, for the target photons of typical en- 
ergy of ~ 10 eV, the major fraction of the absorbed energy 
will be released in the form of ~ 100 GeV electrons. Thus 
secondary pair synchrotron emission must show a spectral 
break in the Suzaku energy band unless one assumes unrea- 
sonably high magnetic fields, B > 1 kG, in the surroundings 
of the gamma-ray emissio n region ( Khangul van et al. f 2008a; 
iBosch-Ramon et al.ll2008ah . 

Figure [5] shows the synchrotron and IC cooling times of 
electrons, as a function of electron energy, calculated for the 
stellar photon density at d = 2 x 10 12 cm and for a magnetic 
field B = 3 G. It can be seen that synchrotron losses domi- 
nate over IC losses at E e > 1 TeV. Note that the TeV gamma- 
ray production takes place in the deep KN regime. This im- 
plies that the cooling time, f coo i = 7/7, of electrons generat- 
ing GeV gamma-rays via IC scattering (Thomson regime) is 
shorter than the cooling time of TeV electrons responsible for 
producing very high-energy gamma-rays (KN). The same ap- 
plies for synchrotron cooling time of multi-TeV electrons that 
produce low-energy (MeV) gamma-rays by synchrotron radi- 
ation. One should therefore expect significantly higher MeV 
(synchrotron) and GeV (IC) fluxes than at keV and TeV en- 
ergies, provided that the acceleration spectrum of electrons 
extends from low energies to very high energies. However, 
in the case of existence of low-energy and very-high-energy 
cutoffs in the acceleration spectrum, the gamma-ray fluxes 
> 10 MeV and at GeV energies would be significantly sup- 
pressed. 

To better understand the energy ranges of the electrons re- 
sponsible for X-ray and gamma-ray production, we show in 
Figure|5]the energy zones of electrons relevant to the Suzaku, 
Fermi, and HESS radiation domains. Note that the recon- 
struction of the average energy of electrons responsible for the 
IC gamma-rays depends only on the well known temperature 
of the companion star T = 3.8 x 10 4 K. The light green zone 
in Figure [5] marked as "Suzaku synchrotron", corresponds to 



electrons responsible for the synchrotron photons produced in 
the energy interval 1 keV < e syn < 40 keV. For a reasonable 
range of magnetic field values, the energy interval of elec- 
trons relevant for Suzaku data overlap on one hand with the 
HESS energy interval, and can overlap with the Fermi one. 
This should allow us, in the case of detection of MeV/GeV 
gamma-rays by Fermi, to considerably reduce the parameter 
space, in particular, to better localize the X- and gamma-ray 
production regions from electromagnetic cascade constraints, 
and derive the broadband energy spectrum of electrons and 
the strength of the magnetic field, both as a function of the 
orbital phase. 

Formally, when X-rays and TeV gamma-rays are produced 
by the same population of very-high-energy electrons, one 
should expect a general correlation between the light curves 
obtained by Suzaku and HESS. In this regard, the similar- 
ity between the Suzaku and HESS light curves seems to be 
natural. However, such an interpretation is not straightfor- 
ward in the sense that two major mechanisms that might 
cause modulation of the TeV gamma-ray signal are related 
to interactions of electrons and gamma-rays with the pho- 
tons of the companion star, i.e. anisotr opic IC scatter- 
ing and photon-photon p air production ( Khangulva n et al. I 
2008a; iDubus et alJ l2008). and thus cannot contribute to the 
X-ray modulation. The X-ray modulation requires periodic 
changes of the strength of the ambient magnetic field or the 
number of relativistic electrons. Note, however, that the 
change of magnetic field would not have a strong impact as 
long as the radiation proceeds in the saturation regime and 
synchrotron losses dominate in the relevant energy interval. 
One would also expect modulation of the synchrotron X-ray 
flux if the energy losses of electrons are dominated by IC scat- 
tering, although in such a case we should observe significantly 
lower X-ray fluxes. 

A more natural reason for the modulation of the syn- 
chrotron fluxes would come from dominantly adiabatic losses. 
The adiabatic cooling of electrons in binary systems can be re- 
alized through complex (magneto)hydrodynamical processes, 
e.g. due to interactions between a black hole jet or a pulsar 
wind with t he dense stellar wind o f a massive companion star 
(see e.g. [Bogovalo v et alJ I2008L iPerucho & Bosch-Ramonl 
2008). The orbital motion could naturally produce the mod- 
ulati on of adiabatic cooling of electrons around the orbit (see 
e.g. Khangulvan et al. 2008b). Note that because of the rel- 
atively small variation of the X-ray flux over the orbit, a fac- 
tor of only two, the requirements for this scenario are quite 
modest. We note that dominant adia batic losses have been in- 
voked by Khangulv an et al. I ( 120071) to explain the variations 
of the X-ray and TeV gamma-ray fluxes from the binary pul- 
sar PSR B 1259-63. 

The detected power-law spectrum of X-rays with photon 
index around T = 1 .5 implies that the established energy spec- 
trum of electrons is also a power-law with index a e ~ 2. 
This agrees well with the hypothesis of dominance of adi- 
abatic losses, because the adiabatic losses do not change 
the initial spectrum of electrons. Thus the required power- 
law index a e ~ 2 implies a reasonable acceleration spectrum 
Q{E e ) oc E~ 2 . Otherwise, in an environment dominated by 
synchrotron losses, the acceleration spectrum should be very 
hard, with a power-law index < 1, or should have an unrea- 
sonably large low-energy cutoff at E > 1 TeV to explain the 
observed X-ray spectra. 

Obviously, adiabatic losses modulate the IC gamma-ray 



Study of LS 5039 with Suzaku 



7 



flux in a similar manner. However, unlike X-rays, the TeV 
gamma-rays suffer si gnificant distor tion due to photon-photon 
absorption (see e.g. Bottcher 2007) and anisotropic IC scat- 
tering with its strong hardening of the gamma-ray spectrum 
dKhangulyan & Aharonianl[2005l) . All this leads to additional 
orbital modulation of the gamma-ray signal, and it is likely 
that these two additional processes are responsible for the 
strong change of gamma-ray flux, much more pronounced 
than that seen in X-rays (see Fig. |2). The Suzaku data pre- 
sented in this paper implies a key additional assumption, 
namely that the accelerated electrons must loose their energy 
adiabatically before they cool radiatively. 

In order to demonstrate that the suggested scenario can 
satisfactorily explain the combined Suzaku X-ray and HESS 
gamma-ray data, we performed calculations of the broad-band 
spectral energy distributions (SEDs) of the synchrotron and 
IC emission, assuming a simple model in which the same 
population of electrons is responsible for both X-rays and 
TeV gamma-rays. We also assumed that the emission re- 
gion has homogeneous physical conditions. This is a reason- 
able assumption given that we deal with very short cooling 
timescales (<C 100 s), thus electrons cannot travel significant 
distances while emitting. 

In the regime dominated by adiabatic energy losses, the 
synchrotron X-ray flux is proportional to f a( j. The X-ray mod- 
ulation seen by Suzaku is then described by the modulation 
of the adiabatic loss rate. In Fig. [6] we show f a d(0) that is 
inferred from the X-ray data. The required adiabatic cool- 
ing timescales are ~ 1 s. Any consistent calculation of the 
adiabatic cooling requires the solution of the corresponding 
hydrodynamical problem, and one needs to know in detail the 
nature of the source. At the present stage, we consider the 
simple example of adiabatic cooling in a relativistically ex- 
panding source. In such a case, the adiabatic loss rate can be 
written as: 7 a d(</>,7) = 7 /fad with f a d ^R/c ~ 3Ru sec, where 
R\[ = R/(10 u cm) is the characteristic size of the source. 
The required variation of the adiabatic cooling is thus re- 
duced to the modulation of the size of the radiation region 
(Rn ~ 0.3-1). The size in turn depends on the external pres- 
sure exerted by, e.g., the stellar wind from the massive star. 
The expected weaker external pressure around apastron im- 
plicitly assumed in our model would be broadly consistent 
with the radial dependence of the wind pressure. 

In Fig. [7] we show the SEDs averaged over the INFC (</> = 
0.45-0.9) and SUPC (4> < 0.45, 4> > 0.9) phase intervals. The 
correspond ing gamma-ray data ha ve been previously reported 
by HESS (lAharonian et al.1 120061) . Since both the absolute 
flux and the energy spectrum of TeV gamma-rays vary rapidly 
with phase, in order to compare the theoretical predictions 
with observations, we should use smaller phase bins, ideally 
speaking with the time intervals Af < 100 s corresponding 
to the characteristic cooling timescales of electrons. Because 
of the lack of the relevant gamma-ray data available to us, 
we here use the X-ray and gamma-ray data integrated over 
A(f> = 0.45. While this compromise does not allow us to per- 
form quantitative studies, it can be used to make a qualitative 
comparison of the model calculations with observational data. 

The theoretical calculations of the SEDs are in a reasonable 
agreement with the observed spectra, though they do not per- 
fectly match the gamma-ray fluxes. One can improve the fits 
by introducing slight phase-dependent changes in the spectra 
of accelerated electrons, but it is beyond the scope of this pa- 
per given the caveat mentioned above. We should also note 




VHE photon index 




Orbital phase 

FIG. 6. — Lightcurves. (a) the phase dependence of adiabatic losses derived 
from the Suzaku data; (b) solid line: calculated 1 TeV gamma-ray fluxes 
together with HESS data points, dotted line: the normalizations at 1 TeV of 
power-law fits for the calculated spectra; (c) calculated photon indices around 
1 TeV together with power-law indices reported by HESS; (d) calculated 1— 
10 keV X-ray fluxes with the reported Suzaku data; and (e) predicted 1-100 
GeV integrated IC flux (note that the peak in the spectrum is at 10 GeV; 
see Fig. 7). The calculations were performed for the following parameters: 
inclination of the orbit = 30 deg; elevation over the orbital plane Zrj = 1.8 X 
10 12 cm; magnetic field strength B = 2.5 G. A power-law injection spectrum 
of electrons Q(E) oc E; 2 (lOGeV < E c < lOOTeV) was assumed. The total 
power injected in relativistic electrons was fixed at the level of 10 37 erg/s. 
The cooled electron spectrum is formed due to radiative (IC and synchrotron) 
and non-radiative (adiabatic) losses, derived from the Suzaku data. 

that the calculations of low energy gamma-rays (in the Fermi 
domain) are performed assuming that the injection spectrum 
of electrons continues down to 10 GeV. If this is not the case, 
the gamma-ray fluxes in the Fermi energy domain may be sig- 
nificantly suppressed. On the other hand, the detection of 
gamma-rays by Fermi would allow us to recover the spec- 
trum of electrons in a very broad energy interval, and thus 
distinguish between different acceleration models. Another 
important feature in this scenario is that a hard synchrotron 
spectrum extending up to a few MeV, is required by the ro- 
bust detection of ~ 10 TeV gamma-rays from the system. A 
future detection of the emission at MeV energies may bring 
important information on the presence of highest energy par- 
ticles in the system. 

The reproduction, at least qualitatively, of the observed 
spectral and temporal features of the nonthermal radiation 
with the simple toy model supports the production of X-ray 
and TeV gamma-rays by the same population of parent par- 
ticles and allows us to derive several principal conclusions. 
The electron energy distribution should be a power-law with 
an almost constant index of a e = 2Tx — 1 ~ 2 to explain the 
X-ray spectra. Note that in an isotropic photon gas when 
the Compton scattering takes place in the deep KN regime 
such an electron distribution results in a quite steep TeV spec- 
trum with photon index T 7 ~ 1 + a e ~ 3. This does not agree 
with HESS observations. However, the anisotropic IC pro- 
vides a remarkable hardening o f the gamma-ray spectrum 
(Khangulvan & Aharonian 2005), in particular, T 7 ~ 2 would 
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be expected for the INFC, and it h as indeed been observed 
using HESS dAharonian et al.ll2006h . We also note that, to 
explain the VHE spectrum at SUPC, we have to assume that 
the emission region is located at a distance «2x 10 12 cm 
from the compact object, in the direction perpendicular to 
the orbital plane. In the standard pulsar scenario, the pro- 
duction region cannot be located far away from the compact 
object, and e ven invoking electromagnetic cascadi ng (see e.g. 
Fig. 16 in iTorres & Sierpowska-Bartosikl 120081) . one can- 
not reproduce the reported fluxes around orbital pha se 0.0 
(iKhangulyan et al. 1 120081 iBosch-Ramon et al.1 l2008al) . We 
have found that the magnetic field strength cannot deviate 
much from a few G. We can also derive a constraint on the 
size of the emission region, imposing a maximum expan- 
sion speed of ~c (the speed of light), and the Hillas criterion 
(Hillas;ll984j), in which the minimum size of a source capa- 
ble of accelerating particles to a given energy E e is R = Rl 
(where R^ = £ e v /300Z?g cm), the Larmor radius. This esti- 
mate yields a size of 10 10 -10 n cm, which agrees quite well 
with the estimate based on the required timescale of adiabatic 
cooling. Note that the requirement of fast adiabatic losses 
imposes a strong constraint on the acceleration rate of elec- 
trons. Indeed, the acceleration timescale can be expressed 
as: f acc = r]R L /c ~ 0.1?7(£e/lTeV)(fi/T G)" 1 s, where -q > 1 
parametrizes the acceleration efficiency. In extreme acceler- 
ators with the maximum possible rate allowed by classical 
electrodynamics r\ = 1. The HESS spectrum provides evi- 
dence of electron acceleration well above 10 TeV. Therefore, 
face < fad ~ 1 s is required at E e = 10 TeV, which translates 
into r\ < 3 for B = 3 G. Thus we arrive at the conclusion that 
an extremely efficient acceleration with r\ < 3 should operate 
in a compact region of R ~ 10 1 1 cm. 

Finally, we would like to emphasize that in the scenario 
described here different radiation energy intervals are char- 
acterized by fundamentally different light curves. While the 
synchrotron X-ray modulation is caused by adiabatic losses, 
the light curve in gamma-rays depends critically, in addition, 
on effects related to interactions with the optical photons of 
the companion star. Two of these effects, photon-photon pair 
production and anisotropic Compton scattering are equally 
important for the formation of the light curve of TeV gamma- 
rays. On the other hand, only the effect of anisotropic Comp- 
ton scattering has an impact on the formation of the light curve 
of GeV gamma-rays. The difference of light curves in the X- 
ray and GeV and TeV gamma-ray intervals in this scenario is 
shown in Figure [6] 

6. SUMMARY 

The Suzaku X-ray satellite has observed LS 5039 for the 
first time with imaging capabilities over one and a half or- 
bits. The Suzaku data show strong modulation of the X-ray 
emission at the orbital period of the system and the X-ray 
spectral data are described by a hard power-law up to 70 keV. 
We found the close correlation of the X-ray and TeV gamma- 
ray light curves, which can be interpreted as evidence of pro- 
duction of these two radiation components by the same elec- 
tron population via synchrotron radiation and IC scattering, 
respectively. Whereas there are at least two reasons for the 
formation of a periodic TeV gamma-ray light curve, both re- 
lated to the interaction with photons from the companion star 
(photon-photon absorption of VHE gamma-rays and IC scat- 
tering in an anisotropic photon field), the modulated X-ray 
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FIG. 7. — Model calculations of the non-thermal radiation spectra of 
LS 5039, averaged over SUPC (low state) and INFC (high state) orbital 
phase intervals for <j> < 0.45, <j> > 0.9 (SUPC, Blue); and for <f> = 0.45-0.9 
(INFC, Red). The flux points of Suzaku X-ray spectra (1^10 keV band) are 
reconstructed with the best-fit function with correction for the interstellar ab- 
sorption. The HESS gamma-ray spectra are taken from Aharonian et.(2006). 
The calculations were performed for the same model parameters as in Fig. [6] 



son for the modulation in X-rays seems to be adiabatic losses 
which should dominate over the radiative (synchrotron and 
IC) losses of electrons. We demonstrate that this assumption 
allows us to explain, at least qualitatively, the spectral and 
temporal characteristics of the combined Suzaku and HESS 
data. In particular, the introduction of adiabatic losses not 
only provides a natural explanation for the rather stable pho- 
ton index of the X-ray spectrum T^- — 1 .5, but also allows one 
to approximately reproduce the TeV gamma-ray spectra. 

The gamma-ray data require a location of the production re- 
gion at the periphery of the binary system at d ~ 10 12 cm. This 
constraint allows a quite robust estimate of the magnetic field 
of a few Gauss to be derived directly from the X/TeV flux ra- 
tio, and an adiabatic loss time of a few seconds to provide the 
dominance of adiabatic losses. In the case of a relativistically 
expanding source, the size of the production region should not 
exceed 10 11 cm. The adiabatic cooling cannot be shorter than 
several seconds, and correspondingly, the size of the produc- 
tion region cannot be much smaller than 10 11 cm, since oth- 
erwise the electrons could not be accelerated up to energies 
beyond 10 TeV, even assuming an extreme acceleration rate 
close to the fundamental limit determined by quantum elec- 
trodynamics. There is little doubt that future simultaneous 
observations of LS 5039 with the Suzaku, Fermi, and HESS 
telescopes will provide key information for understanding the 
nature of this mysterious non-thermal source. 
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